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We present a method of smoothing discrete breakup S-matrix elements calculated by the 
^ ' method of continuum-discretized coupled-channels (CDCC). This smoothing method makes 

it possible to apply CDCC to four-body breakup reactions. The reliability of the smoothing 
method is confirmed for two cases, **Ni(d, pn) at 80 MeV and the El transition of ^He. We 
0^ ' apply CDCC with the smoothing method to a ®He breakup reaction at 22.5 MeV. Multistep 

breakup processes are found to be important. 

§1. Introduction 

Q . Recent developments in radioactive beam experiments have made it possible to 

^ \ study unstable nuclei away from the stability line. Such nuclei have exotic properties, 

^ i ' e.g., the halo structure, ^•'"^^ in which weakly bound valence neutrons extend far from 

a core nucleus. Borromean nuclei such as ^He and ^^Li are typical examples of halo 
\ nuclei and are described well by a three-body (n-|-n-|-core) model. In the scattering 

• of a three-body projectile, it easily breaks up into its constituents, and hence the 

Q^ , reaction should be described as four-body (n-|-n-|-core-|-target) scattering. 

' Thus far, many experiments have been conducted on the scattering of ^He from 

heavy^^"^'^^ to light targets^^"-^^^ at high®-*"^-^) and low incident energies. ^^"'^^'^^^ From 
the theoretical point of view, it is quite difficult to solve four-body scattering exactly. 
00 ! At higher incident energies, therefore, the scattering has been analyzed using ap- 

■ proximate methods such as the Glauber model, ^^-''^^^ adiabatic approximation,^^)' 

multiple scattering expansion^^^ and the four-body DWBA.^^)"^") These methods 
are, however, not applicable to scattering at low incident energies. 
^ ■ One of the most useful and reliable methods for the scattering in a wide range of 

^ ■ incident energies is the method of continuum-discretized coupled channels (CDCC).^^)'^^-* 

In CDCC, the total scattering wave function is expanded in terms of the complete set 
of bound and continuum states of the projectile. The continuum states are classified 
by linear and angular momenta, k and respectively; each of them is truncated at 
a certain value. The fc-continuum is then divided into small bins and the continuum 
states in each bin are averaged into a single state. This procedure of discretization is 
called the average method. The ^-matrix elements calculated with CDCC converge 
as the model space is 

extended. 23). 24) xhe converged CDCC solution is the unper- 
turbed solution of the distorted Faddeev equations, and corrections to the solution 
are negligible within the region of space in which the reaction takes place. 2^)' 

Conventional CDCC based on the average method has been applied only to 
three-body breakup processes in the scattering of two-body projectiles. It is thus 
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called three-body CDCC. For four-body breakup processes in the scattering of three- 
body projectiles, the average method is not feasible, since it requires exact three- 
body continuum states of the projectile that are quite difficult to obtain. This 
problem can be circumvented by the pseudostate discretization method instead of 
the average method. In the pseudostate discretization method, the continuum states 
are replaced with pseudostates obtained by diagonalizing the internal Hamiltonian 
of the projectile in a space spanned by L^-type basis functions. One can adopt the 
Gaussian^"^^'^^^ or transformed harmonic oscillator (THO)^^^ basis as the L^-type 
basis functions. The validity of the pseudostate discretization method was confirmed 
in the scattering of two-body projectile by the good agreement between the CDCC 
solutions obtained by the pseudostate discretization and the average methods. 

The pseudostate discretization method makes CDCC applicable to the scatter- 
ing of three-body projectiles. In fact, four-body CDCC based on the pseudostate 
discretization method with Gaussian^'^^"^^) or THO^^-* basis functions has been suc- 
cessful in describing the elastic scattering of a three-body projectile at not only high 
energies but also low energies near the Coulomb barrier. This shows that the back- 
coupling effects of four-body breakup channels on the elastic scattering, i.e., virtual 
four-body breakup processes, are described well by four-body CDCC. Further devel- 
opment is, however, necessary to apply four-body CDCC to real four-body breakup 
processes, as shown below. 

Let us consider the four-body reaction system shown in Fig. 1. The three-body 
projectile B, which consists of a, b, and c, is incident on the target nucleus A. In 
the final stage of the breakup reaction considered, four particles, a, b, c, and A, are 
emitted. The cross section of the reaction is described by the quintuple differential 
cross section 



where Sa {sb) is the energy of the emitted particle a (b), and Qa, ^b, and 12c are the 
scattering angles of a, b, and c, respectively. This cross section is calculated using the 
true breakup S-matrix elements S{k, K, P) for the transition of B from the ground 
states to the continuum states classified with three momenta; k, K, and P are, 
respectively, conjugate to the coordinates r,y, and R shown in Fig. 1. The single 
differential cross section da/de with respect to the intrinsic energy e of B is obtained 
by integrating Eq. (1-1) over th.6 five Vctrictbles Sa-, Eb, f2a, f^h, and ilc with energy 
and momentum conservations. However, the breakup S-matrix elements S calculated 
using CDCC are discrete functions of e, because of pseudostate discretization, i.e., 
the diagonalization of the internal Hamiltonian of the projectile. Thus, we need a 
smoothing method for deriving continuous breakup S'-matrix elements S{k,K,P) 



In this paper, we propose a new smoothing method for constructing S{k, K, P) 
from S to complete four-body CDCC. First, we test the smoothing method in 
breakup reactions of a two-body projectile and confirm that the breakup S'-matrix 
elements calculated by the smoothing method agree with the "exact" ones, i.e., the 
results with the average method of discretization. Next, we apply the smoothing 
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method to the electric dipole (El) transition of ^He to its continuum states and 
confirm that the calculated El transition strength, with an additional smearing pro- 
cedure concerning the experimental resolution, converges as the model space is ex- 
tended. Finally, we apply the smoothing method to four-body breakup reactions of 
^He by ^''^Bi at 22.5 MeV and show that multistep breakup processes are significant. 

This paper is constructed as follows. In §2, we recapitulate four-body CDCC 
and present a new smoothing method. Test calculations of the smoothing method 
are shown in §3 and its application to the ^He+^^'^Bi scattering at 22.5 MeV is shown 
in §4. Section 5 gives the summary. 

§2. Formulation 

2.1. Four-body CDCC 

We start with the four-body Hamiltonian of the four-body system shown in 
Fig. 1: 

H = Tr + U + Hb, (2-1) 
U = Ua{Ra) + V^°^\Ra) + U.iRb) + + UdRc) + ^''""'(i^c), (2-2) 

HB=Tr+Ty + Vbc + Vca + Kfe, (2-3) 

where R is the relative coordinate between the center of masses (cm) of B and A; we 
assume A to be a structureless and inert nucleus. Note that back-coupling effects of 
inelastic channels concerning the excitation of A are taken into account with the use 
of the optical potential between A and each constituent of the projectile, as will be 
mentioned later. The coordinate R^ (x = a, b, and c) denotes the relative coordinate 
between x and A, and = R, r, and y) is the kinetic-energy operator associated 
with ^. The interaction V-x^' is the potential between x and x', and t/x is the optical 
potential between x and A; the Coulomb part V-^°'^^ is also accurately treated in 
order to describe Coulomb breakup processes. 

The basic assumption of four-body CDCC is that the four-body reaction takes 
place in the model space^^-'''^^^ 
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where i denotes a set of quantum numbers, i.e., the energy index of the pseudostates 
n, the total spin of the projectile /, and its projection on the z-axis m. The pseu- 
dostate ^i(= ^nim) satisfies 

{^i\HB\^,,) = 6u>i^, (2-5) 

where ii is the eigenenergy of ^j. The validity of this assumption is justified by the 
fact that the calculated elastic and total breakup cross sections of the four-body 
scattering converge as the model space is extended.^^^'^^^''^'^^"^^'' One may thus 
regard {^i} as a complete set in describing the reaction process considered. We 
henceforth call {^i} the approximate complete set in this meaning. 

Following the discussion above, four-body CDCC starts with the four-body 
Schrodinger equation in the model space V: 

V[H - E,,,]V\^^^^^) = 0, (2-6) 

where E'tot is the total energy of the system. The four-body wave function tp'^DCC 
is expanded by the approximate complete set {^i}: 

l^'^°^^)=El^-^^)' (2-7) 

i 

where \^i,Xi) = ® |Xi)) and i = denotes the elastic channel and others {i ^ 0) 
the breakup channels. The expansion coefficient |xi) describes the relative motion 
between B (in state and A. The intrinsic energy ii of B and the relative momen- 
tum Pi between B and A satisfy the energy conservation 

Etot = -Ef™ -I- eo = '^^i + (2'^) 

where fi is the reduced mass between B and A, and Ef^ is the incident energy of B 
in the cm system, i.e., E^^ = h^pQ/{2ii). 

Multiplying Eq. (2-6) by from the left leads to a set of coupled differential 
equations for called the CDCC equation, 

[Tr + Ui,i - Ei]\x^) = -Yl U^AX^'), (2-9) 

where Ei = ffP^/{2fi) = Stot — ^i- The coupling potential f/j^j' is defined by 

Ui,i> = {S,\U\Si,). (2-10) 
The CDCC equation (2-9) is solved under the usual boundary condition for {Il\xi) = 

2.2. Pseudostate Discretization Method 

In the pseudostate discretization method,^'^^'^^^'^'^^"^^^ the pseudostates {^i} are 
obtained by diagonalizing the internal Hamiltonian of B in a space spanned 
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by L^-type Gaussian basis functions. The eigenstates of with negative and 
positive energies correspond to the bound state(s) and pseudostates, respectively, 
and the latter are regarded as discretized continuum states as mentioned above. 
Now we consider ^He as a projectile. The eigenstate = ^n./m is then written as 

3 

{rM'^nim) = E^il(r-c,2/c), (2-11) 

c=l 

where c denotes a set of Jacobi coordinates defined in Fig. 2. Each V'ri/m ^ °^ 
the Gaussian basis functions: 

V'ilCr-cJ/c) = j;4'^)"Vi'^)(r-c,2/c), (2-12) 

a 

^^^\r,, y,) = MuM^xrie-^-'h^ce-''^''" 

X "[yK^?r.J®n(%)]^® hi;^^®??i;^^]J^^, (2-13) 

where a is the abbreviation of {i, j, ^, A, yl, 5}; I (A) is the orbital angular momentum 
associated with the coordinate (yc) and rji/2 is the spin wave function of each 
valence neutron, ni or n2- The normalization coefficients Ma and Mj\ are determined 
so as to satisfy (A/'j^rf e"''*'"^ jAAj^r^ e"'''''^) = 1 and {Mjwl e~^^y^^\J\fj\yc e^^^^^ = 1. 
The quantum numbers ^, A, and A are truncated by the upper limit values ^maxi 
Amax, and ylmax; respectively. The total spin S is either or 1. The Gaussian range 
parameters are given geometric progression: 

Ui = \/rl n = n , (2-14) 

A, = l/y|, y:, = yx{yr...lyif-'^^l^'---^y (2-15) 

The states ^nim are antisymmetric under the interchange between n\ and n2, and 
hence they must satisfy ^^^"^ = (— and for c = 3 (— )'*'^'^ = 1. Meanwhile, 
the exchange between each valence neutron and each nucleon in ^He is treated ap- 
proximately by the orthogonality condition model. ^^-^ The eigenenergies e^i of ^He 
and the corresponding expansion coefficients J^a^^ are determined by diagonalizing 
^g_36),37) diagonalization procedure one needs the coordinate transformation 

between the rearrangement channels. Details of this transformation are described in 
Appendix A. 



6 



T. Egami, T. Matsumoto, K. Ogata and M. Yahiro 



2.3. Smoothing Method 

The exact breakup T-matrix elements to the continuum state ip{k, K) of B are 
given by 

rE^(fc, K, P) = K), P\Um, (2-16) 

with \ip{k, K), P) = \ip{k,K)) (S> \P), where k, K, and P are the momenta in 
the asymptotic region associated with the coordinates r, y, and R, respectively; 
\'ip{k,K)) is the exact three-body wave function of B with the energy e satisfying 

[HB-e]mk,K)) =0, (2-17) 

and |P) is the plane wave function satisfying 

[Tn-{Etot-e)]\P) =0. (2-18) 

The exact four-body (a-fb-|-c+A) wave function can be replaced with the corre- 
sponding CDCC wave function ip-CDCC -^[^i^ good accuracy. Inserting the approx- 
imate complete set V of Eq. (2-4) between the bra vector and the operator U on 
the right-hand side of Eq. (2T6), we can obtain the approximate smooth T-matrix 
elements T{k, K, P): 

Tik, K, P) = K)\S,){Si, P,|C/|y^CDCC^ 

i 

= Y,Hk,K)f,, (2-19) 

•i 

where = \^i) \Pi), ^i{k,K) is the smoothing factor defined by 

^,(fc,K) = (V'(fc,K)|<^,), (2-20) 

and Ti is the breakup T-matrix element of CDCC defined by 

T; = (^„A|[/|ip-CDCC^^ (2-21) 

Since the breakup T-matrix elements are proportional to the breakup S-matrix ele- 
ments, Eq. (2-19) is reduced to 

S{k,K,P) = ^J^i{k,K)S^. (2-22) 

i 

The smoothing factor J-'i{k, K) is the overlap between the pseudostate and the 
exact continuum state 'ip{k,K). In principle, the exact continuum state is obtained 
by solving the three-body Schrodinger equation (2-17) under the initial condition 
that the three particles a, b, and c, are incident particles. In practice, however, it is 
quite difficult to do so. Therefore, we approximately solve the equation in the model 
space V of Eq. (2-4): 

= Vli^o) + V —i— — -pypiv^), (2-23) 

e — Ir— ly + l€ 
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where V = Vat + Vhc + Vca and i/jq describes free propagation of the three incident 
particles satisfying 

[Tr + Ty-£] iVo) = 0. (2-24) 

The solution to Eq. (2-23) should converge as the model space is extended. This is 
tested in §3. Multiplying Eq. (2-23) by leads to 

Ti{k, K) = ^i{k, K)+Y^ GijVjkMk, K), (2-25) 

jk 

where 

^i{k,K) = {S,\^Jo) = {i^om, (2-26) 

Vjk = {'P,\V\Sk)- (2-28) 

Since Eq. (2-25) is a set of linear equations for J^i(k,K), one can easily obtain a 
solution for J-'i{k, K) once ^i{k, K), Gij, and Vjk are evaluated. We emphasize here 
that what we need in the present smoothing method is not the exact three-body con- 
tinuum wave function il){k,K) itself but the smoothing factors J^i(k,K). In other 
words, as shown in Eq. (2-20), we need only the overlap between ip{k,K) and 
the latter is an accurate three-body wave function in the model space that is signif- 
icant for the four-body reaction under consideration. Thus, henceforth, we call this 
procedure based on model space truncation the model space smoothing method. The 
six-fold integral over r and y in Eq. (2-26) can be made analytically, and the six-fold 
integral in Eq. (2-27) is reduced to a single integral. These properties markedly sim- 
plify numerical calculations. We show in Appendix B the explicit form of 
i.e., the Fourier transform of and in Appendix C that of Gij. 

We remark that in J^i{k,K), or equivalently, in Eq. (2-23), the rearrangement 
channels of the three particles are fully taken into account. On the other hand, in 
the calculation of scattering processes by four-body CDCC, we neglect four-body 
rearrangement channels, in which at least one constituent of the projectile is bound 
in the target nucleus. The contribution of the rearrangement channels is, however, 
theoretically shown to be negligible in forward angle scattering. ^^^'^^^ Recently, it 
has also been confirmed numerically by comparing the results of three-body CDCC 
for deuteron elastic and breakup processes by ^^C at 56 MeV with the exact solution 
of the Faddeev equation. ^^-^ Such a comparison for four-body scattering processes 
will be very interesting. 

§3. Test Calculations 

3.1. Breakup Reaction of Two-body Projectile 

In this subsection, we consider the ^^Ni{d,pn) reaction at 80 MeV to test the 
validity of the model space smoothing method. In this reaction, the projectile (d) 




Fig. 3. Smoothing factor as a function of k for n = 16. The left (right) panel corresponds to £ = 
{£ — 2). The dots and solid curve denote the results of the exact calculation and those obtained 
by the model space smoothing method, respectively. 



consists of a proton (p) and a neutron (n), and the total system is composed of p, 
n, and A. The Coulomb breakup processes are neglected, i.e., Coulomb potential is 
treated as a function of the relative coordinate R between the cm of d and ^^Ni. 
We adopt the model space smoothing method described in §2.3 to obtain continuous 
breakup 5-matrix elements for this reaction. Note that, in this case, the coordi- 
nate y and its conjugate momentum K do not appear and r denotes the relative 
coordinate between p and n. As for the interaction between p and n, we cite that 
in Ref. 27). The ground and breakup states of the p+n subsystem, classified with 
the linear and angular momenta k and i, respectively, are obtained by diagonalizing 
the Hamiltonian of a subsystem with complex-range Gaussian basis functions. As 
for the parameter set of the basis functions, we take {2A'max = 40, ri = 1.0 fm, 
^max = 20 fm}, where 2iVmax is the total number of the basis functions, and ri and 
rmax are the Gaussian range parameters given in Eq. (2-14). The maximum momen- 
tum /cmax is taken to be 1.3 fm~^, and the number of states is 17 for each of £ = 
and 2. The optical potentials of n-^^Ni and p-^^Ni systems are cited from Ref. 27). 
In solving the CDCC equation, the scattering wave functions are connected with the 
asymptotic forms at Umax — 

30 fm. 

When a two-body projectile such as d is considered, it is possible to obtain an 
exact smoothing factor without making the model space truncation in Eq. (2-23). In 
Fig. 3, the approximate smoothing factor with the model space smoothing method 
(solid curve) is compared with the exact smoothing factor (dots) in the case of n = 16 
with i = (left panel) and 2 (right panel). In each panel, the two results agree very 
well with each other. In fact, the difference between the two is below 1% level. 

Figure 4 shows the squared magnitudes of the breakup S'-matrix elements for 
the total angular momentum J = 17; panels (a), (b), (c), and (d) correspond to 
{i,L) = (0,17), (2,15), (2,17), and (2,19), respectively. Again, the results obtained 
by the model space smoothing method (solid curves) agree very well with the results 
of the exact calculation (dots). 
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Fig. 4. Squared magnitudes of the smoothed breakup S-matrix elements [^(fc)^ as a function of k 
at the total angular momentum J = 17 for the d+''*Ni scattering at 80 MeV. Panels (a), (b), 
(c), and (d) correspond to {l,L) = (0, 17), (2, 15), (2, 17), and (2, 19), respectively. The results 
of the exact calculation (model space smoothing method) are shown by the dots (solid curve). 

3.2. Electric Dipole Transition of ^He to 1" Continuum States 

In the calculation of the ^He+n+n three-body wave functions, we take the same 
Hamiltonian as that in Ref. 31), which is based on the orthogonality condition 
model (OCM).^^^ In the present calculation, however, the OCM potential to ex- 
clude the forbidden state is introduced into only the 0+ state. Consequently, one 
can obtain the 0"*" ground state of ^He while excluding the forbidden state in the 
diagonalization of H^. On the other hand, for the 1~ state, not only the physi- 
cal pseudostates to be regarded as discretized continuum states of the ^He+n+n 
three-body system, but also unphysical eigenstates that have large overlaps with the 
forbidden state are obtained, because of the absence of the OCM potential. We 
exclude the latter by hand, following the variation-before-projection method. This 
rather special treatment of the 1~ state is due to the following fact. In the present 
calculation, as one can see below, so many Gaussian basis functions are used for 
the 1~ state to cover a wide range of internal coordinates concerned. In this case, 
when an OCM potential that has a core of about 10^ MeV appears, the numerical 
calculation of the diagonalization of the ^He internal Hamiltonian becomes unstable. 
Note that the 0"*" pseudostates and 2"*" state of ^He, which are included in Ref. 31), 
are not relevant to the present calculation of the El strength distribution of ^He. 
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Table 1. Parameters of the basis functions for the state of He. 



c 


£ 


A 


A 


S 


^max 


ri 
(fm) 


(fm) 


J max 


J/i 

(fm) 


^max 

(fm) 


3 

1,2 
3 

1,2 
1,2 




1 
1 
1 




1 
1 
1 




1 

1 





1 



1 


10 
10 
10 
10 
10 


0.1 
0.5 
0.1 
0.5 
0.5 


20.0 
20.0 
20.0 
20.0 
20.0 


10 
10 
10 
10 

10 


0.5 
0.5 
0.5 
0.5 
0.5 


20.0 
20.0 
20.0 
20.0 
20.0 


Table II. 
Set I 


Three parameter 


sets (I, II, and III) of the basis funct: 


ions for the 1 


^ state of 


«He. 


c 


£ 


A 


A 


S 


^max 


ri 
(fm) 


^max 

(fm) 


J max 


yi 

(fm) 


^max 

(fm) 


3 

1,2 
1,2 




1 


1 
1 



1 
1 
1 







17 
17 
17 


0.5 
0.5 
0.5 


50.0 
50.0 
50.0 


17 
17 
17 


0.5 
0.5 
0.5 


50.0 
50.0 
50.0 


Set II 


3 

1,2 

1,2 





1 


1 
1 



1 
1 
1 







18 
18 
18 


0.5 
0.5 
0.5 


70.0 
70.0 
70.0 


18 
18 
18 


0.5 
0.5 
0.5 


70.0 
70.0 
70.0 


Set III 


3 

1,2 
1,2 





1 


1 
1 



1 
1 
1 







19 
19 
19 


0.5 
0.5 
0.5 


80.0 
80.0 
80.0 


19 
19 
19 


0.5 
0.5 
0.5 


80.0 
80.0 
80.0 



As for the basis functions taken in the diagonalization of Hb, we adopt the 
real-range Gaussian basis functions. The parameter set for the 0"*" state is shown in 
Table I. For the 1~ state, we take the three parameter sets shown in Table II and see 
the dependence of the El strength distribution on them; the model space described 
by the 1~ pseudostates of ®He taken in the calculation dictates the convergence of 
the El distribution. Note that for the 1~ states, rmax and ymax are taken to be 
50-80 fm to cover the large model space concerned. The S = 1 components of the 
1~ states are found to be negligible and not included in the present calculation. 

The electric dipole (El) strength from the 0^ ground state to the nth 1~ 
pseudostate is 

5(El;n) = iV |(^„i„|Oi^(El)|^g^)p, (3-1) 

where Oi^(El) = r^Yi^iilra) and N = {N„NjAf. The mass number and neutron 
number of the core nucleus are denoted A and Nc, respectively, and Nn is the number 
of valence neutrons; for ^He, A = Q and Nn = Nc = 2. 

Since the final states {'Pnim} of the El transition have discrete energies {in}, 
B(El;n) is a discrete function of the energy e of the ^He+n+n system in the final 
state, as shown in Fig. 5 (the bars); panels (a), (b), and (c) show the results with 
Sets I, II, and III, respectively. B(El;n) should be smoothed in some way to be 
compared with the corresponding experimental data.^^'^'^^ 
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Fig. 5. El strength distribution of ®He as a function of e. In each panel, the bars show _B(El;n) 
from the 0^ ground state to the pseudostates of ®He. Panels (a), (b), and (c) correspond to the 
results of Sets 1, II, and 111, respectively. The solid, dashed, and dot-dashed curves respectively 
show the smeared El strength distributions with F = 0.5, 0.2, and Q.liin + 0.975). 



A possible simple way is to assume a Lorentzian function for the e dependence 
of each S(El;n):39) 

dBCEl) r 1 , , , , 
^ — L = S^- Bm-n). (3-2 

In Fig. 5, the solid (dashed) curves show the results of Eq. (3-2) with a width F of 0.5 
(0.2) MeV. With F = 0.2 MeV, the El strength has an unnatural oscillation, which 
is independent of the parameter set of the basis functions. Even if F increases to 
0.5 MeV, the distribution becomes smooth only in the small e region of e < 3 MeV. 
Also shown by the dot-dashed curves are the results with the energy-dependent width 
F = 0.1 (fn + 0.975) cited from Ref. 9). Again, an unnatural oscillation is found in 
each panel. In all the cases, unphysically, the smeared distribution is finite at e = 0. 
Thus, this simple way of smoothing is not helpful for obtaining a result that can be 
compared with the experimental data.^^'^*^^ 

The exact El strength B^^{E\) is defined by 



(El) = Ar // dkdKY,\{'^im{k,K)\Oi^{El)\^ll)\\ (3-3) 
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Fig. 6. El strength distribution of He as a function of e. The left and right panels correspond 
to the nonsmeared and smeared results, respectively. In each panel, the dotted, dashed, and 
solid curves respectively show the results of Sets I, II, and III. Also shown for comparison by 
the dot-dashed curve is the result of the simple smoothing method with the energy-dependent 
i.e., the dot-dashed curve in Fig. 5(c). 



where ipimik, K) is the exact continuum state of 1 with the momenta k and K. 
The El strength distribution function dB^^(El)/de is obtained using 



de 



NjJdkdK5 (e-^-i^) j;|(V;i^(fc,K)|Oi^(El)|^^^)p. 



(3-4) 



Inserting the approximate complete set V of Eq. (2-4) between the bra vector {tpimik, K) 
and the operator Oi^(El) on the right hand side of Eq. (3-4) leads to the El strength 
distribution dB{El) / de based on the model space smoothing method: 



2 



(3-5) 

where 

A^(El;n) = (^„i„|Oi^(El)|#g^), (3-6) 

Tnl^{k, K) = ii^lmik, K)\Snlm). (3-7) 

Thus, the El strength distribution function dB{Fjl)/de is obtained from the smooth- 
ing factor J-nimik, K) and the discrete El transition matrix element A^(E1; n). Note 
that Al(El;n) is related to the discrete El strength i?(El;n) through Eq. (3-1). 
The main task to obtain dB{Ejl) / de is, therefore, to evaluate the smoothing factor 

nlm{kj 

Figure 6(a) shows the dependence of dB{El)/de on the model space, i.e., the 
three parameter sets in Table II. The dotted, dashed, and solid curves represent the 
results calculated with Sets I, II, and III, respectively. The El strength distribution 
function almost converges as the model space is extended. Strictly speaking, however, 
the convergence is not perfect. This is because the largest model space of the ^He 
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wave function that we can treat in the present calculation, i.e., Set III, is stih not 
sufficient to give a perfect convergence. Note that we here consider the Coulomb 
dipole excitation of ^He. Since the corresponding interaction has a long range, the 
model space of ^He, in the 1~ state in particular, requires quite large values of rmax 
and ymax of the Gaussian basis functions, which considerably increases the number of 
the basis functions. On the other hand, in the nuclear breakup process of two-body 
deuteron shown in §3.1, one can prepare an approximate complete set of deuteron 
quite easily, because not only this breakup potential concerned is short-ranged but 
also deuteron has only one internal coordinate. 

Thus, at this stage the model space smoothing method cannot give a pure the- 
oretical result fully converged owing to computational limits. In order to compare 
the theoretical result with experimental data, however, one may use an additional 
smearing procedure that takes account of experimental resolution. This procedure 
is expressed by 

As for the weight function w{e, e'), we use the Lorentzian distribution with an energy- 
dependent width r = 0.1(e' + 0.975) cited from Ref. 9). Figure 6(b) represents 
the model space dependence of the El strength distribution function smeared with 
Eq. (3-8). A clear convergence is seen in the figure and the result converged is 
consistent with those of the other existing methods. '^'^^^^^^ This is also the case 
if we adopt the Gaussian distribution for tt;(e,e'); the difference in the result from 
that with the Lorentzian distribution is within the thickness of the curve shown in 
Fig. 6(b), at least for e < 10 MeV. Also, as shown in Fig. 6(b), the dot-dashed 
curve represents the result of the simple smoothing method of Eq. (3-2) with Set 
III and the energy-dependent F, i.e., the dot-dashed curve in Fig. 5(c). One sees 
the dot-dashed curve deviates considerably from the solid curve. Thus, the simple 
smoothing method is not accurate as mentioned above. This is found to be also the 
case when the above-mentioned w{£,e') is taken in the simple smoothing method 
using Eq. (3-2). 

Through the tests in §3.1 and 3.2, one can say that the model space smoothing 
method is reliable in the description of experimental data of the breakup reactions 
concerned. Therefore, four-body CDCC based on the pseudostate discretization 
method and the model space smoothing method is expected to be applicable to four- 
body Coulomb and nuclear breakup reactions such as the (^He,^Henn) reaction. 

§4. Application to Four-Body Breakup Cross Section 

In this section, we apply the model space smoothing method to the breakup 
reaction of ^He by ^''^Bi at 22.5 MeV. The four-body Hamiltonian H taken here is 
the same as that in Ref. 32) except that the OCM potential is introduced into only 
the 0"*" state of ^He. Moreover, since this is the first trial of four-body CDCC to 
real four-body breakup reactions, for simplicity, the model space in the present four- 
body CDCC calculation is taken to be smaller than that in Ref. 32); we include only 
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the 0"*" ground state and 1~ pseudostates of ^He in the coupled-channel calculation. 
These states are obtained as in §3.2; the parameter set I is used for the 1~ state. In 
CDCC calculation, the maximum internal energy ^max of ^He is taken to be 7 MeV, 
which results in 170 1~ pseudostates. The scattering waves between ^He and ^'^^Bi 
are numerically integrated up to iimax = 200 and connected with the asymptotic 
forms. The maximum total angular momentum Jmax is taken to be 200. 

The angle-integrated breakup cross section to the continuous breakup channel 
with the momenta k and K is given by 

aik, K) = Y: ^^^^^ I E Sidk, K)\\ (4-1) 

J,L Im 

where Pq is the initial momentum of ^He in the cm system, J is the total angular 
momentum of the four-body system, and Sj^{k, K) is the breakup S'-matrix element 
for the transition from the initial channel to the breakup channel with (/, L, J). Here, 
/ is the total spin of the projectile, and L is the orbital angular momentum associated 
with R. 

Inserting the partial-wave expansion form of Eq. (2-22) into Eq. (4T), one can 

find 

a{k,K) = E ''%2 \Y.^nUk,K)S'^ji^ \ (4-2) 



p , 

J,L ^ nIm 



where 5"^/^ is the discrete breakup S-matrix element for the transition from the initial 
channel to the rath discrete breakup channel with the quantum numbers (I, L, J). 
Thus, the angle-integrated energy spectrum of the breakup reaction is obtained by 

Since we do not obtain a fully converged theoretical result, we calculate, as in §3.2, 
the smeared result with 

^ ' de w{e,e ) — , (4-4) 



de / J de 
, 0.1(e' + 0.975) 1 ^, 

"^"'^) = (.-.')^ + {0.1(e' + 0.975)}2- ^^-'^ 

Figure 7 shows the result of {da/de) as a function of e. The solid curve is the re- 
sult of four-body CDCC and the dashed curve is the result of the one-step calculation 
in which only one-step transition from the elastic channel to the individual breakup 
channels is taken into account. The latter is significantly larger than the former. 
This clearly shows the significant reduction in the size of the breakup cross section 
due to multistep transition. Furthermore, the two calculations give quite different 
e dependences for e > 2 MeV. Therefore, the assumption of a one-step transition is 
not valid at low incident energies. The one-step assumption is often adopted when 
El strength distribution is derived from a measured breakup cross section. Even 
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300 



four-body CDCC 
one-step calculation 




2 3 4 5 6 
f [MeV] 

Fig. 7. Angle-integrated energy spectrum (da/de) of ^EIe+^'"'Bi scattering at 22.5 MeV. The solid 
curve shows the result of the four-body CDCC calculation and the dotted curve shows the result 
of the one-step calculation. These results have been smeared by the same procedure described 
in S3.2. 



though this has been conducted at higher incident energies,^)' ^'^^ where multistep 
processes are expected to be less important, it will be very interesting and impor- 
tant to quantitatively evaluate the contributions of multistep breakup processes in 
this energy region. 

§5. Summary 

We propose a new smoothing method based on model space truncation, i.e., 
the model space smoothing method, for constructing a smooth breakup cross section 
from a discrete one obtained CDCC. This method allows the evaluation of not only 
three-body breakup cross sections but also four-body ones, which is an important 
advantage to the smoothing method established so far. 

The reliability of the model space smoothing method is confirmed in the two 
cases: ^^Ni(d, pn) at 80 MeV and the El transition of ^He. In the former, the new 
model space smoothing method is found to give the same breakup S-matrix elements 
as the "exact" ones based on the previous smoothing method in Ref. 27). In the 
latter, the El strength distribution obtained by the model space smoothing method, 
with the additional smearing procedure of Eq. (3-8), turns out to converge well as 
the model space is extended. 

We apply the model space smoothing method to the ^He scattering on ^'^^Bi at 
22.5 MeV, i.e., near the Coulomb barrier energy, and calculate the angle-integrated 
energy spectrum of the four-body breakup reaction. Comparing the result of the 
four-body CDCC calculation with that of the one-step calculation, we find that the 
one-step calculation is not accurate at low incident energies. 

Thus, the description of four-body breakup reactions with four-body CDCC is 
accomplished by the model space smoothing method proposed in this work. In a 
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forthcoming paper, we will analyze the experimental data on the angle-integrated 
energy spectrum da/de^^'^^^ and the energy and angular correlations'^^ among ^He 
and the two neutrons in the breakup of ^He at lower and intermediate incident 
energies. 
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Appendix A 

Coordinate Transformation Between Rearrangement Channels 



The coordinate transformation from {rb,yi)) to {rc,yc) is defined by 

('^^\ = ( \ ( '^''\ 

W) \ibc ^'bcj \yc) • 

With this transformation, the basis function 'p^alifbiUb) is rewritten as 



(A-1) 



^ I c iJc 



[Y,^{Qr.)(^Yy^{Qy^)]^^ ® [rit/^®r,^^/i]s^\^^ , (A-2) 

where ab is the abbreviation of {ih-,3b-,^b-,\-,Ah^Sh} for channel h and 

r]bc = ^i^llc + Aj,7'L (A-3) 

Cbc = y^Ac + Aj,5'L, (A-4) 

ibc = ^ijbclbc + >^jb^'bcl'bc- (A-5) 

The coefficient (^f,A;,ylfe|£^A^Tylfe)^^^ in Eq. (A-2) is that of the Raynal transformation 
defined by 

J2 (4AfeAK;A',r.lfc),^,r^''+^''-^yn^^;,(^'-J ® YyS^yc)]A, 



(24 + 1)(2A, + 1) £ £ J (2A)!j24-2A)!V mI^^' 



A=0 A=0 
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'4 -A \b-A l[ 
X lt''it~^^LS'iY.{ ^ ^ K) (4 - AOAft - .io|/,o)(ao.io|a;,o) 

4 Aft Au 



where T = \ + A. 



(A-6) 



Appendix B 

Fourier Transform of 'Pnlmi'r',y) 



The Fourier transform ^nlm{k, K) of ^nlmii',y) has the form 

1 



(27r)2 



(27r)3 Z^^" V^a l^c^c 



(B-l) 



with 



vr 



3/2 



vr 



3/2 



2ui. 



ic 



2A,. 



Im 



(B-2) 



Performing the coordinate transformation of Eq. (A-1) leads to 



TT 



3/2 



A/'i.^AG.A, ( — ) ( — ) ( ^ 



vr 



'J6 



3/2 



2A,, 



As 



X ky^'-^K^ 



2 J5i, 



7m 



with 



Vbc 



1 



4m 



1 



4A,, 



be 



where 



7bc = Icb, Sbc = I'cb: I'bc = ^cb, S' 



be = 0^. 



(B-3) 

(B-4) 
(B-5) 
(B-6) 

(B-7) 
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Appendix C 

Matrix Elements of Gij 



In Gij the six-fold integration is reduced to a single one: 



~ X] X] X] ^a}''^ ^a}^^Gaaat,, 
a=l 6=1 aa oil, 

where Oc = {ic,jc,^c, Ac, Ac, Sc} for channel c (=1-3) and Gaaa^ is given by 



SS' 



(2vr)f 



vr vr vr vr 

^ia ^ja ^ib ^3b 



3/2 



2Vi. 



la 



2A 



'3a 



2A 



(C-1) 



X (^;,040|kO)(A;,OA60|«;0) 



xH) 
with 



dre' 



7r(2m + 2K + l)!! 

2m+«;+n+4 



fc=0 
n-m //-/\2fc+K 



(2/c + 2K + 2n + 1)!! /^m 



(2A; + 2k + 1)! 



{{C' + irily){v' + iTfir)-{m 



(C-2) 



/2r = k,,/{2firb), 

V' = r]ab + yi'^^ib), 
2n = £a + h + Xa-T 



jly = h'Ki/i2^iy^), 

C' = C6 + 1/(4A,-J, 
2m = Xb + T — K, 



— Cafe) 



and i = \/2i + 1 and W{a, b, c, d; e, /) is the Racah coefficient. The use of Eqs. (B-3) 
and (C-2) greatly simplifies numerical calculations. 
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